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Abstract 

The bending angle of light is a central quantity in the theory of gravitational lensing. We 
develop an analytical perturbation framework for calculating the bending angle of light rays lensed 
by a Schwarzschild black hole. Using a perturbation parameter given in terms of the gravitational 
radius of the black hole and the light ray's impact parameter, we determine an invariant series for 
the strong-deflection bending angle that extends beyond the standard logarithmic deflection term 
used in the literature. In the process, we discovered an improvement to the standard logarithmic 
deflection term. Our perturbation framework is also used to derive as a consistency check, the 
recently found weak deflection bending angle series. We also reformulate the latter series in terms 
of a more natural invariant perturbation parameter, one that smoothly transitions between the 
weak and strong deflection series. We then compare our invariant strong deflection bending-angle 
series with the numerically integrated exact formal bending angle expression, and find less than 1% 
discrepancy for light rays as far out as twice the critical impact parameter. The paper concludes by 
showing that the strong and weak deflection bending angle series together provide an approximation 
that is within 1% of the exact bending angle value for light rays traversing anywhere between the 
photon sphere and inflnity. 

Keywords: gravitational lensing, black holes 
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I. INTRODUCTION 



One of the early striking predictions of General Relativity is that the weak-deflection 
bending angle of star light grazing the Sun is of the form 

2" 



«Einstem('"o) 



m. 



+ O 



m. 



(First-Order Weak-Deflection) (1) 



which at leading order is twice the value given by Newtonian gravity. Eddington's 1919 
confirmation of the leading term in ([1]) was the first observation of gravitational lensing 
and brought Einstein's new gravitational theory instant scientific acclaim. Gravitational 
lensing in the weak-defiection limit has since been studied extensively, yielding numerous 
applications in astrophysics and cosmology (e.g., Schneider et al. 1992 [2], Fetters et al . 2001 
1], Kochanek et al. 2005 [9]). In addition, over the past two years weak-defiection lensing 
has been employed to create tests accessible to current or near-future instruments, of gravity 
theories such as PPN models and 5-dimensional, string-theory inspired, braneworld gravity 
(Keeton and Fetters 2005-2006 Q,^^). In Q, the first-order weak deflection formula ^ 
was also extended to all orders in m,/ro and re-expressed as an invariant perturbation series 
(since tq is a coordinate dependent quantity [7|). Interested readers may also find a recent 
analysis in of the weak deflection limit. 

In recent years, however, the exciting promise of planned space-borne black hole imaging 
instruments has ignited research activity in the analytical study of lensing in the strong- 
deflection regime (e.g, Virbhadra and Ellis 2000 [9|, Frittelli, Kling, and Newman 2000 lOj, 
Eiroa, Romero, and Torres 2003 (llj. Fetters 2003 (jj], Ferlick 2004 (isj, Bozza, Capozziello, 
De Luca, lovane, Mancini, Scarpetta, and Sereno 2001-2005 14|). 

For a Schwarzschild black hole of physical mass M, the spacetime geometry in the vicinity 
of the photon sphere at radius r = 3m,, where m, = GM/c^ is the black hole's gravitational 
radius, is revealed through the resulting strong-deflection gravitational lensing. In 1959, 
Darwin [l^ computed the first-order term of the bending angle of light traversing deep 
inside the black hole's potential — i.e., close to the photon sphere: 

36(2 - V3) m. 



"Da 



(ro) 



-Ti + 2 loe 



+ O [h'] , (First-Order Strong-Deflection) (2) 



ro - 3m, 

where h' = 1 — (3m,)/ro. He also showed analytically that near the photon sphere there 
are two families of relativistic images, which are images determined by light rays that loop 



around the photon sphere at least once before reaching the observer. Other authors have 



confirmed t 
minet 1979 



lis strong-deflection multi-looping lensing effect (e.g., Atkinson 1965 |16[|, Lu 



17|, Chandrasekhar 1982 |l8|, Ohanian 1987 [l9|, several recent authors j9t|-jl4|). 



These studies were based on evaluating the lowest-order term (out from the photon sphere) 
of the light ray's strong-deflection bending angle. Equation ([2]) is the well-known leading 
logarithmic deflection term. 

In this paper, we develop a perturbative framework that allows us to generalize Darwin's 
strong deflection result ([2]) to any order in h' . Surprisingly, we also found that the leading 
oga rithmic deflection term employed in the literature can be improved. Earlier studies (e.g.. 



15l . p. 188], 18|, p. 132]) arrived at the leading logarithmic expression by a perturbation 



scheme that seems to combine higher and lower order terms. (We leave a definite assessment 
to the judgment of the reader — see the Appendix.) By re-doing the perturbation theory 
and being careful to compare only terms of the same order, we obtain an improvement 
(i.e., more accurate expression) to the leading logarithmic deflection term. Furthermore, 
since ro is coordinate dependent Tj, we re-formulate our strong-deflection bending angle in 
terms of a coordinate-independent series. In particular, we compute this invariant series 
explicitly to 3rd-order in the perturbation parameter b' = 1 — bc/b, where b^. = 3v^m, is 
the critical impact parameter. Our perturbation framework was also used to compute the 
weak-deflection bending-angle series directly and we found it to be in complete agreement 
with the expansion found recently in ^. Finally, we show that our invariant bending-angle 
series is in excellent agreement with the numerically computed exact formal expression for 
the bending angle. This is done for both the strong and weak deflection limits, and the span 
from the photon sphere to infinity. These results should be applicable to analytical lensing 
studies across these regimes and serve as a limiting case to check bending angle results in 
spacetime geometries generalizing the Schwarzschild metric. 

The outline of the paper is as follows: Section [IT] expresses the light's bending angle in 
a formal exact expression involving a difference of elliptic integrals of the first kind. In 
Section IIIIl we expand the strong-deflection bending angle in terms of an invariant series 
going outward from the the photon sphere. This section includes the improvement to the 
logarithmic term. Finally, Section HVl gives a numerical comparison between the perturbative 
and exact bending angles across the range from the photon sphere to infinity. 
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II. FORMAL EXACT STRONG-DEFLECTION BENDING ANGLE 



A Schwarzschild black hole is the unique static, spherically symmetric, asymptotically flat 
vacuum solution of the Einstein equation. The metric is given in Schwarzschild coordinates 

{t,r,9,(P) by 

ds' = - (l - ^) + (l + ^) ' + f2 {de' + sin^ 9 d<p'), (3) 

where t = ct and m, = GM/c^ (gravitational radius) with r physical time and M the 
physical mass of the black hole at the origin. 

Consider a standard gravitational lensing situation where a point source and observer lie 
in the asymptotically fiat region. In a typical lensing scenario, the source and observer are on 
opposite sides of the black hole. However, this restriction can be lifted in strong-deflection 
lensing. Suppose that the source is close to the optical axis passing through the observer 
and black hole. By spherical symmetry, it suffices to choose the source-to-observer light rays 
as lying in the equatorial plane = vr/2). The Euler-Lagrange equations yield that the 
light rays are governed by (e.g., [4]): 

^ 1 

(4) 



drj - (1 - 2m,/r)/r^ 

where b = \L/E\ is the impact parameter with L and E the respective angular momentum 
and energy invariants of the light ray. Setting u = 1/r, re-write as 



du^ ^ 



— (l-2m.M) 



u 



= 2my~u^ + ^. (5) 

This cubic polynomial has a maximum of two positive roots and at most one negative root. 
Writing ([5]) as 

B{u) = 2m, (m — ui){u — U2){u — u^), 

we consider the case of one negative root ui and two distinct positive roots U2 and u^. The 
three roots, given in terms of an intermediate constant Q that allows us to line up the roots 
in the order Ui < U2 < U3 are given by (e.g., p. 130 [isl): 

tq - 2m, - Q 1 ro - 2m, + Q 

Ml = , U2 = — , M3 = . 

4m,ro ro 4m,ro 



Here ro is the light ray's distance of closest approach, which is determined from (e.g., Eq. 
(12) of fl): 

^3 

(6) 



ro - 2m. 

By comparing the coefficients in B{u) to those in the original polynomial in equation ([5]), 
we obtain the following two relations between Q and the quantities b, m,,ro: 



which is equivalent to 



8m, fg 



g2 = (ro-2m.)(ro + 6m.). 



n 



The bending angle of the lensed light ray is given by (e.g., Eq. (20) of [^]): 



a 



-y/2m,(-u - ui){u - U2){u - Ms) 



— TT 




2 [•^/''■Q 



du 



TT 



Jo v^(u - Ui){u2 - U){U3 - U) 

Split the above integral into two parts to make the lower limit equal to the smallest root ui: 



a 




m. 



U2 



du 



du 



n. (7) 



/«i a/(m - Ui){u2 - U){U3 - u) Jui \/{u- Ml) (^2 - u) (% - u) _ 

Now, the integrals in ([7]) can be realized as elliptic integrals of the first kind (see Byrd 
and Friedman 21| for an introduction to elliptic integrals): 



a 




m. 



2F(^i,A;) 2F(^2,A;) 



. V^S - Ui ^/U3 - Ui 

where -^(^'j, k) is an incomplete elliptic integral of the first kind with amplitudes 



and modulus 



Explicitly, 



T ^ T • -1 / 

^1 = -, ^2 = sm \ , 

2 \ U2 — Ui 



2 U2-U1 



U3 - Ui 



. _i /Q + 2m. -ro 2 Q-ro + Qm, 

\&2 = sm \l— z , k - 



Q + 6m, - ro' 



2Q 
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Hence, the exact bending angle simplifies to 

a = A^[Kik) - Fi^.k)] - TT, (8) 

where K{k) and -F(\E', k) are the complete and incomplete elliptic integrals of the first kind, 
respectively, and = \E'2- As a check of ([8]), note that in the limit when m, 0, we obtain 
Q = vq, k = and \1/ = 7r/4, which in turn imply K{k) = tt/2 and k) = tt/4 to give us 
zero deflection as expected. 

It is important to add that the exact bending angle expression ([8]) serves mainly as a 
formal expression. The latter has to be evaluated to obtain explicit analytical and physical 
properties about the nature of the bending angle. Equations ([1]) due to Einstein and ([2]) 
due to Darwin are the first-order evaluations of ([8]) in the weak and strong deflection limits, 
respectively. The challenge of course is in evaluating ([8]) beyond those terms. The weak- 
defiection series in terms of the impact parameter out to many orders beyond ([I]) was 
found recently in Q]. We shall now determine the strong-deflection series beyond ([2]), re- 
derive the weak series result in {4], and reformulate the weak series in terms of a new 
perturbation parameter to allow a seamless comparison covering the span from the strong 
to weak deflection limits — i.e., from the photon sphere to infinity. 



III. EXPANSION OF BENDING ANGLE BEYOND THE PHOTON SPHERE 



The photon sphere is defined by the radius r = 3m,, which marks an unstable photon 
orbit. Light rays that cross within the photon sphere are captured by the black hole (e.g., 
0, 3)- Using the relations given in Section [III we can see that exactly on the photon 
sphere the impact parameter invariant is given by the critical value be = 3^/3m,. Photon 
orbits of interest to us are between tq = 3m, and tq = 00, and indeed our focus will be on 
the region closer to the photon sphere. 

We first express h' in terms of b'. Equation ([6]) is a cubic in r^ that is readily solved to 
yield: 

/ — 3-\/3 m. 



2b 

7! 



cos 



COS 



(9) 



3 \^ b 

The quantities tq and b of course have very different physical meaning for the light rays. 
Relative to an inertial observer at infinity, the quantity tq is the distance of closest approach 
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to the center of the black hole, while the impact parameter b is the perpendicular distance 
from the black hole's center to the asymptotic tangent line to the light ray converging at 
the observer. Overall, the quantity tq approaches b as we extend into regions well beyond 
the photon sphere, but near the photon sphere, the values of vq and b are different. 

To seamlessly traverse from regions near the photon to those at infinity, a natural choice 
of invariant parameter is 

,/ bf, 

which ranges from at the photon sphere to 1 at infinity (asymptotically fiat region). 
A. AfRne Perturbative Form for Bending Angle 

Our goal is to show that the strong-deflection bending angle can be expressed as an 
"affine perturbation" series in b'. More precisely, define an affine perturbation series about 
a function g as 



where Ai and Bi are constants with p and q positive rational numbers. We shall demonstrate 
in Section [ill DI that the bending angle has an invariant affine perturbation series of the form 



where Aq, and pi are numerical constants. Note that ffTTj) is not a Taylor series expan- 
sion because of the appearance of the logarithmic term. However, we shall see that this 
logarithmic term is not exactly ([2]). 

B. Bending Angle Series Beyond the Logarithmic Term 

We now consider the exact bending angle in the region around the photon sphere by 
expanding out from the photon sphere using h = m./rg. The expression for the bending 
angle can now be rewritten as a function of h using the following useful relations: 



fix) = (Ao + . . . + ApxP + ■■■) g{x) + (5o + ■ ■ ■ + B^x'^ + ■■■), 



(10) 




(11) 



Q = roV(l-2/i)(l + 6/i) 
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v/(l - 2h){l + Qh) - (1 - Qh) 
2y/{l - 2h){l + 6h) 



sm 



v/(l -2/i)(l + 6/i) - (1 - 2h) 



v/(l -2/i)(l + 6/i) - (1 - 6/i) ' 
As tq increases from 3m, to oo, the parameter h goes from 1/3 to 0, which corresponds 
to the invariant impact parameter increasing from 3 ^/Sm, to oo. Intuitively, when the 
bending angle is computed in the regime near b = 3-\/3m,, we speak of strong-deflection 
(since we shall show in Section [III Bl that the bending angle can become arbitrarily large), 
while weak-deflection will refer to regions with large b. 

Before proceeding with the substitution of the above quantities in the expression for 
bending angle, we present an example of the standard series expansions for a complete 
elliptic integral of the first kind (e.g., see page 298 in 2l|]) in terms of the modulus k: 



K{k) 



71 

2 



l + + -k' + —k' + 
4 64 256 



or, in terms of the complementary modulus k' = \/\ — k'^ 



where fon = and br 



m=0 
b ' 2 



-1/2 
m 



J2m 



(12) 



(13) 



^m— 1 



+ 



These, along with other similar expansions for 



2m(2m-l) ' 

elliptic integrals, are used in the following: After re-expressing all quantities in terms of h, 
we use the above-mentioned series expansions for elliptic integrals to obtain an expression 
for the bending angle ([H]) in powers of h. Expanding around h = (region at infinity) yields 



a(h) =4:h + 
+ 



-4 



15 
T 



-TT 



122 15 



7783 3465 



2 

21397 



TT + 



3465 , , , 
-130 + ^—TT 1 h'^ 



64 



310695 



-TT 



10 16 J ' \ 6 ■ 256 +^[(^)]- (14) 

This is in exact agreement with the weak- deflection bending angle series found in [4]. 

To obtain the strong-deflection bending angle, we expand a{h) out from the photon 
sphere h = 1/3. It is then convenient to work in terms of h' = 1 — 3h. The coefficient 



^/tq/Q of the bending angle (fT5l) is given in terms of h' by 



3 + 4(1 -/i')/i' 



1/4 



while the modulus and amplitude are 



,2 V3 - ^VSh' + J3 + A(l - h')h' 
k = — 



2^y3 + 4:{1 - h')h' 



and 



Since the quantity 



= sin~^ 



3 + 6h' - V3^/3 + 4(1 - h')h' 
12h' ' 



1 - k 



2 -VS + 2V3h' + ^3 + 4(1 - h')h' 



2V3 + 4(l-/i')/i' 

approaches zero as /i' — 0, it is convenient to re- write ([H]) in terms of k': 



a 



K{Vl - k'^) - F(^, Vl - A;'2) 



7r. 



(15) 



An /I'-series expansion for the elliptic integrals in (fT5|) can then be obtained by first expanding 
in terms of fc'^ and then expanding in h' . 

Remark: It is important to point out that if you want Mathematica to compute 



the functions K{-\/l — k'"^) and -F(\E', yl — k^) as given in our notation, then the 
commands are EllipticK[l — A;'^] and EllipticF[\E', 1 — fc'^], respectively. It is incor- 
rect to use EllipticK[\/l - and EllipticF[^, Vl - k''^] because Mathematica defines 
EllipticF[^, m] = (1 - m sin^ 0)-i/2 and EllipticK[m] = EllipticF[7r/2, m]. 
The bending angle series for f|T5|) in terms of h' is then 



d + TT 



8 tanh"^(V2 - Vs) + log(144) - 21og(/i') 



+ 



2h' 



+ - 



+ 



'2 



1 



162 L 
1 

'5184L 
1 



h 



/3 



38880 



-35 + 6 v/3 + 60 tanh"^ ( V2 - Vs) + 15 log(12) - 15 log(/i') 
295 - 36 - 360 tanh"^(y2 - v^) - 90 log(12) + 90 log(/i') 
-19637 + 3420 + 27720 tanh'^( ^2 - v^) + 6930 log(12) - 6930 \og{h') 
211679 - 34200 - 277200 tanh"^(V2 - Vs) - 69300 log(12) 



,'5 



+ 



1 



559872 L 



- 69300 \og{h') 
-5580415 + 989964 v^S + 7456680 tanh"^(V2 - v^) 

+ 1864170 log(12) - 1864170 \og{h') \ h'^ + O [{h'y] 
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The series ( !T6|) can be simplified significantly using the following: 



8 tanh-^(v/2 - ^3) +log(144) - 21og(/i') = 2 log 



12(2-^3) 



60 tanh-^(y2 - ^3) + 151og(12) - 151og(/i') = 151og 



360 tanh"^(v^ - Vs) + 90 log(12) - 90 \og{h') = 90 log 



12(2 - v^) 



12(2 - V3) 



27720 tanh"^(V2 - V3) + 69301og(12) - 69301og(/i') = 6930 log 



12(2 - V3) 



where we made use of tanh"^(V2 - Vs) = j log(2 - ^3). The bending angle series can 
now be expressed more simply as follows: 

' 12(2-^3) 



^, , 12 2- V3 \ 2h' 1 
a = -vr + 21og , + — + — 

^ ' h' / 3 18 



+ 



+ 



1 



162 
1 

5184 



295 - 36 - 90 log 



'l2(2- V3)^ 



h' 



6^3 -35 + 15 log 



h' 



f\2 



■19637 + 3420^3 + 6930 log ( ^^^^ , 

\ h' 



{h'Y + o [{h'f] . 



(17) 



We can extend this series arbitrarily far, but the expressions become cumbersome and there 
is no need to do so for our later invariant analysis. Observe that a becomes arbitrarily large 
as /i' — > due to the first logarithmic term (and since the other logarithms are dominated 
by the given powers of h'). This is the reason for the terminology strong deflection for light 
rays passing near the photon sphere. 



C. Comparison with Darwin's Logarithmic Bending Angle 



The lowest order /I'-term in ( IT71) is given by 
a + TT , /'l2(2-V3) 



log 



h' 



+ ^ih'] = logl^^^^^l + 0[h']. (18) 
tq — om. 



At first glance the reader may think that this term is the well-known logarithmic term found 
by Darwin 1959 15(] and employed commonly in the literature (e.g., Eqs. (262) and (268) 
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in 18|, p. 132]). However, Darwin's result is actually 



2 



+ 7r 



log 



36(2- y3)m. 



r-Q - 3m, 



) 



(19) 



Equations (|T8|) and (fT9|) are not identical! If the quantity vq in the numerator of equation 
(ITS]) were replaced by 3m,, we obtain Darwin's result, but this is not a legitimate substitution 
since the denominator would become zero. However, in the limit where tq approaches the 
photon sphere's radius 3m,, equations (ITS]) and (IT^ would both become infinitely large and 
so will be close in value in that limit. 

In Section [TV] we shall show explicitly how ([18]) is a better approximation than ([19]) . 
For the convenience of the reader, we review in detail in the Appendix how (119]) is derived. 
Note that in the derivation of ( IT8l) . we were careful throughout to compare only terms at 
the same order, which allowed us to read off the lowest order term directly from the series 
expansion. 

D. Invariant Bending Angle to Third Order Beyond Logarithmic Term 

The strong-deflection bending angle series (IT7]) is coordinate dependent since it is given 
in terms of the distance tq of closest approach. We now determine the bending angle in 
terms of the invariant quantity 




where be (critical impact parameter) is given by be = sVSm,. The quantity ranges over [0, 1] 
since b increases outwards from the photon sphere at 6 = 6c to infinity. 

In the Appendix, we derive the Darwin term (IT9]) and show that it is equivalent to the 
following — see equation (\A9\} : 



To express out bending angle in terms of the invariant b' , we first write h' in terms of b' 
using the relationship to get 



+ vr = log 



216(7-4^3) jl-b') 
b' 



(20) 




(21) 
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Series expanding (12T]) in terms of h' gives 



h' 



+ 3. ib'r 



7 



54^6 



A3/2 



+ 



243 



'A 2 



91 



3888^6 



(b 



/^5/2 



2717 



ib' 



1 7/2 



{b'Y + O [{h'f'^] . (22) 



6561 ' 419904v^ " ' ' 59049 

For the weak deflection bending angle series ([HD, note that it is in a non- invariant form. 
To obtain an invariant expression, convert to the parameter h' by using h = {1 — h')/3 and 
then employing the series for h' in ( !22l) . This yields 

128 



a{b') 



^(1 - b') + —(1 - bf + ^1 

3v^^ ^36^ ^ 243v^ 



1-6') 



+ - b'f + - b'f + o[{i- yf] 



5184 " ' ' 10935^3 
The first term is the Einstein term ([T]) given in terms of b' by 

4(1 - b') 



(23) 



Q^Einstcin(^') 



3^3 



(24) 



Turning to our strong-deflection bending angle expansion, insert the series (1221) into the 
/I'-series ( |T71) for strong deflection. This will yield a series in terms of the invariant b': 

~216{7-4^/3)' 

a{b') = — TT + log 



216(7-4^3) 



1^ 



b' 



-879 + 236 V3 + 205 log 



+ 



216 (7-4v^) 



1296 

321590 + 90588 ^/3 + 68145 log 



ib 



i\2 



216(7-4v^) 
b' 



The lowest order term 



629856 



^lowest order + 71" = log 



ib 



/\3 



O [ib'f] . (25) 



216 (7 - 4^3) 



(26) 



is similar, but not identical, to the Darwin logarithmic term (]20l) since (1261) takes into account 
the improvement to Darwin's result. The other terms in the series have not appeared in the 
literature before. Although the series expansion of h' in fl22|) involves fractional powers of 
6', there are no fractional power terms in the bending angle. 
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As promised in Section IIII At we rewrite (l25l) as an affine perturbation series about the 
logarithm: 

a{b') = (ao + a,h' + a,{h'f + a,{h'f + ■ • ■ ) log 

+ (po + Pih' + P2{h'f + p,{h'f + ■■■), (27) 
where Aq = 216(7 — 4v^) and the a's and p's are constants given by 

— TT 

-17 + 4v^ 
18 

-879 + 236^3 
1296 

-321590 + 90588^3 
629856 ■ 

The terminologies below will be employed for the terms of the affine series ( 1271) : 

Oth-order: ctq log (Ao/&') + po 

Ist-order: (ctq + aih' ) log (Ao/6') + (po + Pi&') 

2nd-order: ((Jq + a^h' + a2{h' f ) log (Ao/6') + (po + Pi&' + P2(&')' ) 

3rd-order: {a^ + fXife' + a^ih' f + ^3(6')') log (Ao/6') + (po + Pi&' + P2{h'f + p3(6')') • 

We have chosen to truncate at 3rd-order because this order will already probe accurately as 
far as to twice the critical impact parameter. These issues are taken up in the next section. 



(To = 1 pQ- 

5 

18 
205 

(72 = P2 ■ 

1296 ^ 
68145 
~ 629856 



IV. COMPARISON OF PERTURBATIVE AND EXACT BENDING ANGLES 

This section will give a numerical comparison of the formal exact bending angle with the 
zeroth, 1st-, 2nd- and 3rd-order affine corrections to the logarithmic term. The comparison 

(^exact — Oi\ 
J X 100%. 
'-^exact / 

In Figured] is plotted the percentage discrepancy between the Darwin term fl20|) written in 
terms of b', our Oth-order term fll8p . and the exact numerical value ([H]), which is represented 
by the horizontal axis. The Oth-order term is closer to the exact value and comes within 1% 
of the exact value at b' ~ 0.03. The deviations increase for larger values of b'. 

If we include terms through to 3rd-order in the strong deflection expression (1271) . then 
unlike the case in Figured], the series comes to within 1% of the exact value for regions well 
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2.50% 




b' 



FIG. 1: Percentage Discrepancy between Darwin's logarithmic term, our Oth-order logarithmic 
term, and the exact numerical value represented by the horizonatal axis. Our Oth-order term is a 
better approximation to the exact value. 

past b' ^ 0.03. Figure [2] shows that the regions could be as far out as to roughly twice the 
critical impact parameter (i.e., b ~ 26c oi V ~ 0.4705) and still be within a 1% discrepancy. 
Beyond this point the discrepancy continues to increase and will not be shown. 

For the 5th-order weak deflection bending angle series (123|) and the Einstein term (12^ . 
Figure [3] shows the percentage discrepancies until fl23|) reaches a 1% deviation. This occurs 
at y 0.4705 and has a larger deviation for greater b' values. 

Figure H] shows both the 3rd-order strong deflection fl27|) and 5th-order weak deflection 
(125]) plotted alongside the exact result. The same comparison, in terms of the percentage 
discrepancy is presented in Figure [5l This graph illustrates that if (1271) is used on the interval 
< 6' < 0.4705 and fl23l) on the interval 0.4705 < 6' < 1, then together these two series 
deviate by at most 1% from the exact value. In addition, the 1% maximum deviation occurs 
at b' ~ 0.4705. These two series can then yield bending angle information rather accurately 
from the photon sphere to infinity. 

V. CONCLUSIONS 

An analytical perturbation framework for calculating the bending angle of light rays 
traversing the gravitational field of a Schwarzschild black hole was given. We expressed 
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FIG. 2: Top graph: Percentage discrepancy between the 0th- to 3rd-order terms, the Darwin 
logarithmic term, and the exact result given by the horizontal axis. The 3rd-order result deviates 
by at most 1% from the exact value with the maximum deviation at b' ~ 0.4705. Bottom graph: 
A close-up of the same graph is plotted to show clearly the higher order corrections in more detail. 



the strong-deflection bending angle in an invariant affine perturbation series in b' about a 
logarithmic function. Our logarithmic expression was shown to be a better approximation 
to the exact bending angle than the logarithmic one found by Darwin, which is commonly 
used in the literature. We also derived the known weak deflection bending angle series as a 
consistency check of our framework. The latter series was reformulated in terms of the more 
natural perturbation parameter b, which smoothly transitions from strong deflection near 
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FIG. 3: Percentage discrepancy corresponding to various orders of the weak deHection bending 
angle. The horizontal axis corresponds to the exact value. The 5th-order expansion has a 1% 
discrepancy at b' « 0.04705. 



2.5 



Exact 

* Strong Deflection (Third Order) 
■ Weak Deflection (Fifth Order) 




FIG. 4: The 3rd-order strong deHection and the 5th-order weak deflection are plotted alongside 
the numerically integrated exact formal bending angle. 



the photon sphere to weak deflection at infinity. Comparison was then given of our invariant 
strong deflection bending-angle series with the numericaUy integrated exact formal bending 
angle expression. We found less than 1% discrepancy for fight rays as far out as twice the 
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FIG. 5: Percentage discrepancy for 3rd-order strong deflection and 5th-order weak deEection. The 
horizontal axis corresponds to the exact value. The two discrepancies criss-cross essentially at 1% 
for b' ^ 0.4705. 



critical impact parameter. This was followed by a further comparison with our invariant 
form of the weak deflection series. It was found that taken together, the two series yield an 
approximation that is within 1% of the exact bending angle value for light rays traversing 
anywhere between the photon sphere and infinity. 
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APPENDIX A: DERIVATION OF THE DARWIN TERM 



For the convenience of the reader, we present here the standard bending angle derivation 



as m 



isl . p. 132]. The formal exact bending angle a is given ([8]) by 

1 / \ 1/2 



-(7r + d)=2(^^) [Kik)-Fi^,k)] 
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where 



Q = A/(ro - 2m,)(ro + 6m,), P 



Q -ro + 6m, 

2Q ' 



\I/ = sin ^ 



/g + 2m, -rp 
Q + 6m, - To' 



For the reader's convenience, the correspondence between our notation and that in 
132] is as follows: (tt + d)/2 = 0oo, a = Q, P = tq, D = b, and Dc = be- 
The perturbation scheme in [iSi, p. 132] is defined by setting 

ro = m,(3 + 5). 
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P- 



(Al) 



Note that 6 does not remain under the value 1 as one moves away from the photon sphere. 
The value of 6 is assumed to be small compared to 3. This is important for the approxima- 
tions to follow. Series expanding and keeping terms at Ist-order in 6 yields 



Q = v/(ro-2m,)(ro + 6m,) = ^/9ml + 105ml + 6^ml = 3m, 
Consequently, 



105 6^' 

1 + 



1/2 



Similarly, 



Q = m, 3 + 



55 



Q -ro + 6m, 

2Q 

m, (3 + f ) - m,(3 + 5) + 6m. 



2m, (3 + f ) 



which gives 



Q-6 + 



9 4 

k'^ = -6. 



56 



56 



-1 



9 



(A2) 



Note that the quantity 6 is now related to the complementary modulus of the elliptic function 
allowing for a series expansion in terms of k'. If we consider the limit as /c' — we get 
p. 132]: 



K{k) = K{Vl-k'') log (^^^ = i log (^p) = ^ log ( 



^)=log6-ilog(5) 



and 



F{^,k) = VI - k' 



log 



Vs + i 
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since \1/ sin ^(l/-\/3)- Taken together, these results yield 



(tt + a) 



log 



6^ V3(V3- 1) 



log 



(A3) 



2(v^+l)2 

In the next step, the 6 is solved for in terms of tq and m, using the first-order scheme (lAip . 
and then substituted in the second-order 5^ in ( 1A3I) to get the leading term in equation (268) 
p. 132]: 

^6^73(73-1)2" 



of 



18 



-(vr + a) 



^log 



log 



2(v^+l)2 

36(2 - 73) m. 



^log 



73 (rp - 3m.)2 

2 



(A4) 



ro - 3m, 

It would seem that the first-order scheme flAip should have been extended to 5^ originally 
if second-order terms in 6 would be considered in the analysis. However, we leave it to 
the judgment of the reader to decide whether the mixing of first and second order terms is 
appropriate in the above derivation. Finally, note that equation (1A4I) is the Darwin term 
quoted at the beginning of the paper — see ([2]). 

We can re-express (IA4I) in terms of m, and the impact parameter b. By equation ([6]), we 
have 

3/2 

b = , (A5) 
Vro - 3m, 

Insert the first-order 5 equation ( lAll) into (1A5I) and expand to second-order in the pertur- 
bation parameter 5: 



m3(3 + 5)3 ^ r- 73 ,2 

'1., o =373m, + ^m,52 
m, (3 -I- o) — 2m, 2 



(A6) 



Writing flA6P to second-order in b yields the result in equation (263) of 

, , 73m,52 

b-br. = . 



(A7) 



Once again, we leave it to the reader to decide whether the mixing of first and second order 
terms is appropriate in the above derivation of (1A7|) . Now, substituting 

ro - 3m, 



5 



m. 



from (jAT]) into (jAT]) yields 



m. 



ro - 3m, 



73 m. 
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Using 



we find from (]A4p tliat 



(v^+1) 



TT + d = lo£ 



36^(2 - Vlf 



m. 



3m, 



log 



648^ (VS- l)^ m. 
(73 + 1)' [h-h,) 



(A8) 



Equations flA4l) and flASP are among the common forms used in the hterature [20(] , and both 
are based on combining first- and second-order terms in 5. Finally, we can also express the 
Darwin term using the variable h' . Since 

3^3 m. _l-h' 

b-bc ~ V ' 

we see that equation (lASIl is equivalent to 

216(7 -4v^) (1-6') 



TT + d = lo£ 



(A9) 



Our study in Section IIII Bl carries out the perturbation analysis for obtaining a consis- 
tently, matching terms of the same order. This yields a different expression for the leading 
logarithmic term, one that is more accurate. 
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